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Abstract. We consider a logarithmically correlated random energy model, namely a model for 
directed polymers on a Cayley tree, which was introduced by Derrida and Spohn. We prove 
^ ' asymptotic properties of a generating function of the partition function of the model by studying 
^ a discrete time analogy of the KPP-equation - thus translating Bramson's work on the KPP- 
CN . equation into a discrete time case. We also discuss connections to extreme value statistics of a 
bjQ' branching random walk and a rescaled multiplicative cascade measure beyond the critical point. 
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1. Introduction 



In [TU] Derrida and Spohn introduced a model for a directed polymer on a disordered Cayley 
tree. They put i.i.d. random potentials on each bond of the tree, considered self avoiding walks 
oj on the tree and the partition function 2{t) = where the sum is over the paths u 

which are self avoiding and of length t, /3 is the inverse temperature and E{oj) is the sum of the 
^ potentials at each bond the walk crosses. They then argued that a suitable generating function 
d for Z, {Gt{x) = E(exp(— e~'^'^^(t)))) satisfies the following non-linear integro-difference equation 

^: r 

■ Gt+,{x)= / p{y)Gt{x + y)''dy, (1.1) 

> . 

^ with initial data Go{x) = exp(— e~^^). Here p is the density of the distribution of the random 
[ potentials and K is the order of the tree: each site apart from the root (which has a single 
^ ■ neighbor) has K + 1 neighbors. 

vn ! 

O ■ The main physical interest in this model is that it is simple enough that one can analyze it in 
^ ', a fair amount of detail, but it is also rich enough to contain non-trivial logarithmic correlations 
• • \ between the energies of the paths and hopefully some universal properties of more complicated dis- 
. ^ ' ordered systems with such correlations. Another way to describe this model is that it is a random 
^ . energy model with logarithmic correlations. As anonther reference to the model, we direct the 
. reader to jl3], where the relationship between random multifractal measures and logarithmically 
correlated random energy models are discussed. In particular, it is explained here why the model 
is logarithmically correlated. 

Instead of analyzing this model in great depth, Derrida and Spohn conjectured that the system 
behaves similarly to a continuum one, where the walks are replaced by Brownian motion and the 
splitting in the tree happens at a random exponentially distributed time. They then argued that 
for = 2, the corresponding quantity Gt{x) satisfies the KPP-equation 

dtGt{x) = \dlGtix) + Gt{xf - Gt{x). (1.2) 
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This equation has been studied extensively by Bramson [5]. His results imply that there is a 
function m^{t) so that Gt{x + m^{t)) converges to a traveling wave solution gp (i.e., for a certain 
c(/3), gp{x — c{l3)t) is a solution to (11. 2p ). One particular phenomenon emerging from Bramson's 
analysis is that there is a phase transition in the system. At a certain critical temperature the 
system freezes. This can be seen for example from the form of c(/3) and g/3{x) as well as the 
asymptotics of m^{t): these all become independent of /3 for large enough /3. 

This freezing seems to be a phenomenon occurring in a wide range of disordered systems (see 
[11] for a further discussion and references). Indeed some non-rigorous work by Carpentier and 
Le Doussal [8j suggests that this freezing is something occurring quite generally in systems with 
logarithmically correlated disorder. 

The discrete model is also related to many purely mathematical questions. For example, at 
zero temperature, only the lowest energy configuration is relevant. So the /? = oo case is closely 
related to the question of extreme value statistics, i.e., finding the distribution of min{Vi, Vjv} 
in the N ^ oo limit, where the Vi have logarithmic correlations described by the model. In the 
case where there are no correlations, this distribution is well known and it is known to extend 
to a large class of correlated random variables. Our analysis will imply that the correlations in 
our model are beyond this universality class. In |14) . there is an interesting conjecture about the 
exact form of the extreme value statistics of a certain logarithmically correlated system. 

Instead of considering each self avoiding walk on the tree separately, one can consider them to 
be a single branching random walk. This allows one to define some useful martingales. Branching 
random walks have been considered more generally and in great detail by Biggins and others (see 
e.g. f5]). The /3 = oo case and the problem of extreme value statistics can be interpreted as the 
problem of finding the distribution of the minimum of a branching random walk (or in our case 
equivalently the maximum since we shall be dealing with symmetric random variables). This is 
a problem that has been studied recently for quite general branching random walks (see e.g. |7j 
and yj). Moreover, this interpretation of the maximum of a branching random walk can be used 
to study the maximum of the discrete two dimensional Gaussian Free Field (see [4J and [6j). 

Another interesting problem related to the model is that one can use the energies of the paths to 
construct random measures on hypercubes - so called multiplicative cascade measures (see [18j and 
[2] for information about multiplicative cascades). For K = 2", one splits the ra-dimensional unit 
hypercube into equal sized hypercubes at stage s. At stage t one gives a hypercube the weight 
^-PEi-^t -^j^gj^g j^^ jg ^]-^g gnergy of a path one identifies with a sequence of nested hypercubes and 
the term is there so that the expected volume of the unit hypercube is 1. Several things are 
known about such measures ( |18( |2j). First of all, a weak limit (as t — )■ oo) exists almost surely. 
The limit measure has positive total mass if and only if /3 < a/2 \ogK and it has no atoms almost 
surely. An interesting question is if we can modify the measure in some simple way so that a limit 
would exist also in the /3 > a/2 log K case. Moreover, if the limit exists, does it have atoms? These 
questions are closely related to the notion of multifractality discussed in [13]. As [6] suggests a 
relationship with branching random walks and the discrete two dimensional Gaussian Free Field, 
one might suspect that these multiplicative cascade measures (and their modified versions) are 
related to the measures of quantum gravity considered in |11] . 

Our primary goal will be to show that in the case that the potentials are standard Gaussians, 
the discrete case indeed behaves as expected, i.e., there is a function m^{t) so that Gt{x + m^{t)) 
converges uniformly to a function which is a solution to a stationary version of the recursion 
relation f ll.ip . As in the continuum case, the shift needed to obtain a non-trivial limit is given by 
m^{t) = G^^(|). The stationary version of the recursion relation will turn out to be 
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w{x)=l -^e-^y'w{x + y + c{l3))^dy, (1.3) 

'2/T 



where 

\ V21ogK, for /? > V21ogK ^ ■ ' 

We shall also call w a traveling wave, c{j3) its speed and the equation it satisfies the stationary 
equation. 

As noted in [T?], from the point of view of studying universality classes of extreme value statis- 
tics, the asymptotic behavior of w is important. In our case, 

lim €^''{1 - w{x)) = C (1.5) 



for some C > and for /3 > a/2 log K, 

lim le^^^^l - w{x)) = C (1.6) 

for some C > 0. Moreover, we shall show that under certain restrictions of the initial data, the 
solution to the stationary equation is unique up to translations. 

We note that the stationary equation enjoys a certain high-temperature self-duality. In the 
high-temperature regime, the equation is of the form 

^x) = y_^-^e ^yw^x + y + - + ^j dy. 



This equation is clearly invariant under the mapping (3^^ = While there certainly IS no 

duality between the physics of the high- and low-temperature regimes, this curious formal duality is 
suspected to be related to physical properties of the model. Indeed, in [15] it was noticed that this 
type of high-temperature self-duality occurs in some more complicated logarithmically correlated 
random energy models (where verifying may not be quite a s simple) and it was conjectured that 
this duality property is intimately related to the freezing occurring. Some further support for this 
conjecture was found in [16j . 

From the point of view of the convergence of the random measure mentioned above, one of the 
first questions to ask would be, does AtZ{t) converge, where At is some deterministic normaliza- 
tion. Since Gt{x + m^(t)) = E(exp(e-^(^+™''W)Z(t))) converges, the only possible normalization 
would be asymptotically a multiple of e~^"^'^^^\ For this, the precise form of m^{t) may be im- 
portant. This could also be relevant for studying the 2-dimensional discrete Gaussian Free Field 
since as seen in [B], m°°{t) is related to the expectation of the maximum of the two dimensional 
discrete Gaussian Free Field. Our result for the asymptotic form of m^{t) is 
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c( /3)t + C '(l), for 15 < v /21ogK 

^{t) = { V2WK t - log t + 0(1) for (] = VThjK . (1.7) 
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While it is intuitively rather clear that the discrete system should behave as the continuum one 
and in the /? = oo case this has been showed even for more general non-linearities f7[ fT], a written 
argument in the generality we are considering seems to be missing. Reading through Bramson's 
work on the continuum case, one notices that many of the arguments he uses work in the discrete 
case with minor modifications. Indeed for the part on m^{t), we shall not go over all of the 
technical details that would be formally identical to those found in [5j, but we shall reproduce the 
main argument in the discrete time language and provide the additional results that might not be 
immediately obvious. Moreover, we shall prove the convergence of Gt{x + m^{t)) independently 
of knowing the precise asymptotic behavior of m^{t). This along with the rather small class of 
initial data we are interested allows one to cut a few corners in following Bramson's reasoning. 

Our proof for the convergence of Gt{x + m^{t)) will be rather different from Bramson's for 
yJ2 log K < j3 < oo and we shall make use of arguments used in various areas concerning similar 
problems in discrete and continuous time. In fact, a secondary goal of this note is to collect different 
kinds of arguments and references to various areas which seem to have been independently working 
on similar problems with different kinds of approaches. The point of this being that to study more 
difficult problems such as the existence of the limit measure discussed above, the distribution of 
the partition function of the two dimensional Gaussian Free Field or problems related to quantum 
gravity, a wider range of tools could be useful. 

When proving convergence, one of our main tools will be a generalized maximum principle type 
result which is a discrete version of one used by Bramson. Some of our arguments will follow 
Bramson's approach and some follow |20| . where a similar recursion relation is studied, but the 
density of the random variables has compact support. Another important tool we shall need is 
a family of martingales related to the branching random walk (see e.g. |3J for a more general 
discussion of such martingales). Their properties in the framework of branching diffusions have 
been studied in [21J , which was applied to the KPP-equation in [T7] . For the asymptotic behavior 
of the solution of the stationary equation, we shall rely on work by Durrett and Liggett [12]. The 
study of m^{t) follows the work of Bramson closely and our main tool will be a discrete time 
Feynman-Kac formula and the analysis of a discrete time Brownian bridge. 



2. Tools for demonstrating convergence: a branching random walk and a 

maximum principle. 

In this section, we shall go over some basic results related to a branching random walk and 
a generalized maximum principle for a certain class of integral operators. A lot of the results 
related to the branching random walk have been found for more general branching random walks 
by Biggins and others (see e.g. [3]). As already mentioned, our discussion about the branching 
random walk will rely on work in [21j and [17J . 

The branching random walk we are interested in is defined in the following manner. We start 
with a particle located at some position x. This particle takes a random step to x + V , where V 
is normally distributed with zero mean and unit variance. After this step, the particle splits into 
K new particles {K is fixed) all of which are located aX x + V . After the splitting, one unit of 
time has elapsed. Each of these particles then behaves as the initial one and independent of the 
others. So at time t (an integer), we have particles. They are grouped into K^~^ clusters of 
K particles. Let us write Xfc(t), k = 1, iT* for the locations of the particles at time t. The 
indexing is so that Xi{t), ...,Xi^(t) are in the same cluster, X^+iit), ...,X2K{t) are in the same 
cluster and so on. 
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The fundamental objects we shall use are the random variables Zj^^t) = e I^(^k{t)+c{i3)t) ^ 

The self-similar structure of the branching random walk gives a useful decomposition: Zp{t + s) = 
'Y^^^^e'^'^'^''^^^^^''^^*^ z'^l^\s ,t) . Here for each t, zj^\s,t) are independent copies of Zjs^s) and the 
corresponding branching random walks start from the origin (implying Z^^\o,t) = 1). They are 
also independent of the process up to time t. Let us write {J-'t}t for the filtration of the branching 
random walk. 

The following results concerning the branching random walk are either direct calculations or 
their proofs are simple modifications of those found in [21j and [9]. 

Lemma 2.1. For /3 < v^21ogX, Zjs, dpZj^ and d'^Zp are martingales with respect to the 
filtration of the branching random walk. 

Lemma 2.2. For /3 < y^21og~^, Zj^ is uniformly integrable. 

Lemma 2.3. Let L{t) = minfcXfc(t). Then — t- — ■\/2 log K almost surely as t — )■ oo. 
Lemma 2.4. \im.t^oo{L{t) + ^/2\ogKt) = oo almost surely. 

Let us write p{y) = To prove convergence in the /3 > a/2 log K case, we shall need 

the following generalized maximum principle. 

Lemma 2.5. Let G] and G^ be given by the recursion relation f 1 1.1 1) and have initial data 
Gq and Gq (measurable and between zero and one). Let us assume that the initial data has the 
following property: there is a point Xq G M so that Gq{x) > Gq{x) for x > Xq and Gl{x) < Gl{x) 
for X < Xo- Then for all t > 1, there is a point Xt G [— oo, oo] so that G^{x) > Gl{x) for x > Xt 
and Gf{x) < Gj{x) for x < Xt- Moreover, if Xt G {— oo, oo} for some t, then Xg = Xt for all s > t. 

Proof: A simple induction takes care of the case Xt G {—00,00} for some t. We then note 
that due to the form of the recursion relation and the fact that the initial data is measurable and 
bounded, G] are differentiable for t > 1. Let us now assume that there is a point G M so that 
Gl{x) > Gl{x) for X > Xt and Gf{x) < Gl{x) for x < Xt- Continuity implies that either there is 
a point a G M so that Gf_^_i{a) = Gl_^_i{a) or xt+i G {—00, 00} so let us assume that such a point 
a does exist. Our goal is to show that Xt+i = a satisfies the desired conditions. Integrating by 
parts, 

Gl^'ia) - GlJia) = J p{y)y {Gl{a + yf - G\{a + yf) dy. 

If Xt < a, then 

Gl,'{a) - Gl,'{a) > [ p{y)y (G?(a + y)"" - G]{a + yf) dy 

jR\{xt-a,0) 

+ {xt - a) f p{y) {G]{a + y)^ - G\{a + y)^) dy. 

J xt—a 

Gf^i{a) = Gl_^_i{a) implies that 

[' p{y) {Gl{a + yr - G\{a + y)^) dy = - [ p{y) {Gl{a + y)^ - G]{a + y)^) dy. 

Jxt-a jR\{xt-a,0) 
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Thus 



Gl,\a) - Gl,\a) > [ p{y){y - (x^ - a)) (G?(a + yf - Gl{a + yf) dy. 

JM.\{xt-a,0) 

Since Xt < a, the definition of Xf imphes that the integrand is positive in the integration region, 
so G^^^'(a) — Gl_^_i{a) > 0. Similar reasoning shows that this result holds for > a as well. 

Thus by induction, if for any t > there is a point G M so that Gf{xt) = Gj{xt), then Gf — G] is 
strictly increasing at this point. Now if there was a point x > Xt^o that Gf{x) < Gl{x), continuity 
and differentiability would imply that for some x* G {xt,x), Gf{x*) — Gl{x*) = and Gf — G\ 
would be decreasing at this point. According to our reasoning, this is impossible so G^i^x) > Gl{x) 
for X > Xf. 

If there was a point x > Xt so that Gf{x) = Gl{x), Gf — G] would be strictly increasing at this 
point. This means that Gf — G} would be negative in some open interval {x — 6,x), with 6 > 0. 
This is contrary to Gf{x) > Gl{x) for all x > Xf. Thus actually Gf{x) > Gl{x) for all x > x^. 

Since Gf — G] is strictly increasing at xt, Gf{x) — Gl{x) < for x G {xt — S,Xt) for some 6 > 0. 
If there were some x < Xt so that Gf{x) — G'J(x) > 0, our previous reasoning would imply that 
Gf{y) — Gl{y) > for y > x, which is contrary to Gf{y) — G\{y) < for ?/ G (xf — 5, Xj). ^ 

Remark 2.6. Looking at the proof, it is easy to see that if we change > to > and < to <, the 
lemma will still be true. 

This concludes the discussion about our main tools for proving convergence. 



3. Convergence for /3 < a/2 log K 

In this section, we shall first consider the existence and asymptotics of traveling waves in 
the c > a/2 log K case and then demonstrate convergence for the case where the initial data is 
asymptotically 1 — Ge~^^ (e-g-? exp(— e~^^)), with j3 < yJ2 \ogK. Convergence will also imply 
uniqueness of the traveling waves (uniqueness up to a translation). Many of the proofs are discrete 
time versions of those in |17J. We will begin by showing that for each c > a/2 log there exists 
a traveling wave with speed c. 

Lemma 3.1. For each c > yJ2 log iC, there exists a function : M — )■ [0, 1] so that w is 
increasing, w{—oo) = 0, w{oo) = 1 and 

w{x) = J p{y)w{x + y + c)^ dy. 

Proof: According to Lemma 2.2, the positive martingales Zpit) = e"'^*-"^''*-*^^^*-'^-'*-' are 

uniformly integrable for /3 < ^/l log K. Since they are positive, they converge. Let us write 
Zfi[oo) for the limit. We also define 

Mfiit) =E {e-^^^'^'^\7i) . 

bmce is a uniformly integrable martingale converging to e almost 

surely. Let 
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One can show (see j21]) that Z^(oo) G (0, oo) almost surely for /3 < v^21og~^ so Wf3{—oo) = 0, 
w^(oo) = 1 and Wp is increasing. Thus we only need to show that is a traveling wave with 
speed c{/3). Using the decomposition of + s) and passing to the limit, we have Z^(oo) = 

Y^^^-^^e''^^^''^^'''^^''^'^'^'^ Zj^\oo,t). Since the zj^'' are independent and identically distributed, we see 
that 

= J] E (exp (-e-'^(^^W+'=('^)*)Z^(oo)) | J"*) 
if' 

k=l 

So ''^/3(^fc('^) + c{P)t) is a uniformly integrable martingale, which implies that 

w^ix) = E^ (^wpiXkit) + c(/3)f) j 

for each t. Setting t = 1, we have 

Wi3{x) = / p{y)w(3{x + y + c{l3))^dy, 



i.e. is a traveling wave with speed c(/3). Now for each c > a/2 log K, we can find a /3 < ^/2^ogK 
so that c = c(/3). □ 

The asymptotic behavior of a traveling wave in both cases (c > a/2 log K and c = a/ 2 log -ftT 
follows from a result by Durrett and Liggett [l_2j. 

Lemma 3.2. If w is a traveling wave of speed c = c(/3) > y/2 log K, w{—oo) = and w{oo) = 1, 
then lima;_!.oo e^^(l — = C for some C > 0. 

For the proof, see |12j . Theorem 2.18 a). For the reader interested in going over the proof, 
we shall provide a short dictionary of what the different quantities appearing in [12] look like 
in our case. First of all, Wi = W2 = ... = Wk = e~^^^^'^^^^^ for all i, where y is a standard 
Gaussian. Our solution to the stationary equation corresponds to their fixed point of a smoothing 
transformation via w{x) = (p{e~^^). Moreover, in our case the function v{a) is given by v{a) = 
logX — |c(/3)^ + — c(/3))^. Thus we see that for (3 < v^2 log K, we are in the case where 

v'{a) < if v{a) = and for /3 > y/2 log K we are in the case v'{a) = when v{a) = 0. Finally, the 
associated random walks are quite simple in our case. The increments of the walks are Gaussian 
and for (3 < ^1 log X, the increments have positive expectation while for fi > ^1 log X, they are 
centered. 



We note that if is a traveling wave, then ty( 
the previous lemma implies that for each /3 < 
wave of speed c(/3) so that e^^{\ — w{x)) ^ C 



■ + a) is another traveling wave for each a. Thus 
log K and each C" > we can find a traveling 
as X — 7- 00. 



Our proof of convergence follows that of Lui's |2U]. Following Bramson's arguments, a complete 
classification of initial data for which convergence occurs is possible, but it requires more work. 
Before we demonstrate convergence, we need the following lemma. 



Lemma 3.3. Let Go{x) > Go{x) — Ae for some ^4 > and (3 > 0. Then for any c > y/2 log K 

Gn{x + nc) >Gn{x + nc) - ^e-^^'+"^(^('^)-"). 
Proof: The proof is by induction. We have 

Gn+i{x + {n+ l)c) - Gn+i{x +{n + l)c) = J p{x + c-y) (Cniy + ncf - Gn{y + nc)^) dy 

< KAe""^^"^^^-^^ j p{x + c- y)e-^ydy 

Theorem 3.4. Let e'^^{l — Go{x)) — i'C>Oasx— i-oo and let w be traveling wave of speed c(/3) 
satisfying e^^{l — w{x)) — > C as x — > oo and w{—oo) = 0. Then Gn{x + c{(3)n) — )■ w{x) uniformly 
on sets of the form [a, oo) for any a G M. Moreover, if Go is increasing, then the convergence is 
uniform on M. 

Proof: For any 5 > 0, we have 

x^oo 1 — w{x ± 0) x^oo eP'-^^°>[l — w{x ± 0)) 

Since e~^^ < 1 < e^^, we can find a number Ls so that 1 — w{x + S) < 1 — Go{x) < 1 — w{x — 5), 
i.e., w{x — S) < Go{x) < w{x + 5), for x > Ls- For any 6 > 0, set Ab = sup^<2,^ e^^(l — Go{x)) and 
— sup2.<2,^ e^^(l — wlx + 6)). These definitions imply that 

w{x + 5)> Go{x) - Ebc"^^ 

and 



Go{x) > w{x -6)- Ake-^\ 
Using the previous lemma and the fact that if Gq{x) = w{x) then Gn{x + nc) = w{x), we have 

w{x - 5) - yi^e-^^'+"^('=(^)-'=(^)) < Gn{x + nc{P)) < w{x + 6) + B^,e-^''+''^^<^^-<'^^\ 

Let e > 0. Choosing h G \/2 log K) implies c{h) < c{j3). Thus for any a G M, we can take n 
so large that niax{Sbe-''^+"''W^)-'=(^)), A6e-''^+"''(^(^)-^('^)} < f for x > a. So we see that for x > a 
and large enough n, 

\'w{x) — Gn{x + c{l5)n)\ ^ 2 ~'~ '^^^('^(^) ~ "^(^ ~ 5), w(x + 5) — w{x)). 

We note that the equation of w implies that w is smooth and that < 1. Thus w is uniformly 

continuous and we can take 6 so small that meix{w{x) — w{x — 6),w{x + 6) — w{x)) < | and 

\w{x) — Gn{x + c(/?)n)| < e 
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for X > a and large enough n. So we see that Gn{x + c{j3)n) — )■ w{x) uniformly on sets of the 
form [a^oo). If Go is increasing, then G„ is increasing for each n. Let e > and set a so that 
w{x) < I for X < a. We then take n so large that + nc{[5)) — w{x)\ < | for x > a. Thus 

Gn{a + c{j3)n) < e. Since G„ is increasing, Gn{x + nc{l3)) < e foi x < a. Since G„ and w are 
non-negative, we see that |G„(x + nc{l3)) — < e for x < a and we have uniform convergence 

on M.^ 

Corollary 3.5. For each G > there is only one traveling wave of speed c > ^/2 \ogK so that 
w{—oo) = and e^'^(l — w{x)) — ?■ G as a; — ?■ oo. Moreover, this traveling wave is increasing. 

Proof: In the proof of the previous result, we did not fix which traveling wave w we are using 
apart from fixing G. Since any such w would then be the limit of the sequence G„(a; + c{j3)n), 
we see that there can be only one such w. Taking Gq{x) = exp{—e~^^), a simple induction shows 
that G[{x) > for all x so we have a sequence of increasing functions converging to w so the limit 
must be increasing. ^ 

This concludes our discussion about the /3 < ^/2\ogK case. 



4. Convergence for Heaviside initial data 

We will now focus on Heaviside initial data, i.e., the /3 = oo case. This will be important to 
us in the next section when we prove convergence in the /3 > a/2 logK case. To get started, we 
note that it is simple to check that G^^ : (0, 1) — i- R is well defined for t > 1. Thus we can define 
myit) = Gi\y). 

Convergence of Gt{x + my{t)) follows from Lemma 2.5: 

Lemma 4.1. Let Gt be given by the recursion relation flLip with Heaviside initial data and 
let y G (0, 1) be fixed. Then the limit function Wy{x) = lim^^oo Gt{x + my{t)) exists. 

Proof: Let to e N \ {0} be fixed. Let us set Gj(a;) = Gt+i{x + my{tQ + 1)) and Gf{x) = 
Gt{x + my{to)). Now < Gi(2;) < 1 for all x so we see that Gq{x) < Gl{x) for x < — my(to) and 
Gq{x) > Gq(x) for x > — my(to), so the conditions of Lemma 2.5 are met. 

Moreover, 

GliO) = Gi„(m,(to)) = y = G,„+iK(to + 1)) = G,\(0). 

As to is arbitrary. Lemma 2.5 then implies that for x > 0, {Gt{x + my{t)))t is decreasing and 
bounded from below by so there must be a limit which we call Wy{x). In a similar manner, we 
see that {Gt{my{t) + x))t is increasing for x < and constant for x = so the limiting function 
: R — )■ M exists. ^ 

From now on we shall fix y = | and write mi{t) = m{t) as well as w = wi. We note that w is 

increasing and w{fd) = |. We shall also write Am{t) = m{t + 1) — m{t). To prove that the limit w 
is a traveling wave, we need some simple properties of m. For example, linearizing the recursion 
of 1 — Gt, one can check that there is a constant G so that m{t) < V21og Kt - 2''^^^ + G. 
Moreover, a simple argument using the form of the recursion relation and the fact that Gt{x+m(t)) 
increases to w{x) for x < 0, implies that (Am(t))( is bounded from below. This is in fact enough 
to show that w is a. traveling wave for some c. 
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Lemma 4.2. There is a unique c so that the hmit function w is a travehng wave with speed c. 

Proof: We first note that wc can find a subsequence of {Am{t))t that converges to some finite 
value. Otherwise m{t) < log Kt — 2~2 ^\o^k + ^ would be violated, since we know (Am(t))t 
to be bounded from below. Let the limit of this subsequence be c. 

From the recursion relation, one can check that < < 1 for all x and t. Thus we have 

\Gt^j^x{x + fn{tk)) - w{x - c)| < \Am{tk) - c| + 1^*^+1(2; - c + m{tk + 1)) - w{x - c)| 0, 
as A; 00 for each fixed x. On the other hand 

Gt^+i{x + m{tk)) = j p{y)Gt^{x + y + m{tk))^ dy ^ J p{y)w{x + y)^dy 
as A; — >■ 00. We conclude that w satisfies 



w{x -c) = J p{y)w{x + y)^dy, 



which is precisely the equation we wanted. To show that the value c does not depend on the 

subsequence we picked, let us assume that there were two such values c and c'. Since w is 
increasing, Gi{x) > for all x and Gt{x + m{t)) increases to w{x) for x < 0, t«(a;) > for all x. 
Also w\x) = J p{y)Kw{x + y + c)^~^w\x + y + c)dy so for a given x, w'{x) can be zero only if 
w is a constant. If w were a constant function, we would have w{x) = w{Q) = \, which does not 
satisfy the equation for w. Thus w is strictly increasing. Since 



w{x ~ ^) ~ j + y)^dy = w{x — c'). 

and w is strictly increasing, c — d.^ 

It follows from the recursion relation that < x/nl and < Vn! for all n, x 

and t > 1, which implies that all of these functions are entire. Using some basic results from 

complex analysis, it then follows that the convergence to w is uniform on M. Using this uniform 
convergence and the recursion relation, one can argue that w{x — Am{t)) w{x — c), which 
implies that Am(t) — )■ c and c < a/2 log K. 

To show that c = v^2log^, we will make use of the branching random walk. 

Lemma 4.3. W{t) = Y[k=i w{Xk{t) + ct) is a martingale with respect to the branching random 
walk and w{x) = E^{W{t)) for all t. 

Proof: Decomposing the product into the product over the clusters and the product over par- 
ticles inside each cluster and using independence we have 
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E%W{t + l)\Tt) = 11 p{y)w{Xk{t) + ct + y + cfdy 
k=i 

= \[w{Xu{t) + ct) 

k=l 

= W{t). 

Since is a martingale, W{W{t)) = E^(l^(0)) = w{x). ^ 
Lemma 4.4. c = a/2 log K. 

Proof: Let us assume that c < a/2 log K. The martingale W that we introduced in the previous 
lemma is positive so it converges. Let us denote the limit by W{oo). Moreover, it is bounded 
above by one so it is also uniformly integrable and w{x) = W{W{oo)). On the other hand 

AT* 

< W{t) = Y[w{Xk{t) + ct) < w{L{t) + ct). 

k=l 

According to Lemma 2.3, L{t) + ct — ?■ — oo almost surely as t — )■ oo. Since w{x) — )■ as a: — t- — oo, 
this implies that W{oo) = almost surely and w{x) = W{W{oo)) = for all x, which is contrary 
to our knowledge of w being a function increasing from zero to one. Thus c = a/2 log K. ^ 

For the asymptotic behavior of w, we rely on [121 again (Theorem 2.18 b)) as in the case of Lemma 
3.2. 

Lemma 4.5. Let w he a traveling wave with speed c = a/2 log K. Then - — — — _). (j 
as x — 7- oo for some C > 0. 

This along with our discussion about the branching random walk gives the uniqueness of the 
traveling waves. 

Lemma 4.6. Every traveling wave w with speed c = a/2 log K satisfying w{—oo) = and 
w{oo) = 1, is given by a translation of the limit of the Heaviside case w. 

Proof: Let : M — )■ (0, 1) be any non-trivial solution to the equation 

w{x) = J p{y)w{x + y + c)^dy 
with w{—oo) = and w{oo) = 1. Also let 

Wy{t) = l[w{Xk{t) + ct + y). 

k=l 

As in Lemma 4.3, one can show that is a uniformly integrable martingale which converges (to 
say ly^(oo)). By Lemma 4.5 there is a x G M so that w(x) ~ 1 — xe~'^^^^^\ Thus for each fixed y 
we have 
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lim 



xe 



-c{x+y+x) 



x^oo —\ogw{x + y) ' 
i.e. for each fixed y G R and e > we can find a D e M so that ior x > D 

-c{x+y+x) 



1 -e < 



xe 



— \ogw{x + y) 



< 1 + e. 



Using Lemma 2.4, we see that taking t large enough, L{t) + ct > D almost surely. This means 
that these inequalities hold when we set x to be Xk{t) + ct for any k. Then summing over all k 
we obtain 



(1 - e)(- \og{Wy{t))) < e-^(^+*) J2(Xk{t) + ct)e-^(^'=W+ct) 



k=l 



This means that 



< {l + e){-logWy{t)). 



d 



Wy(oo) = exp e-"(^+^) lim TT7T^/3(t) 

' t-^oo Op 



Since is uniformly integrable, we have 



d 



w{y) = E°(V^^(oo)) = E° exp e-'^^^z+s) Xmi ^ZB{t) 

' ' t-^-oo Op 



/9=c, 



Hence 



- x) = E" ( exp ( e-''" \\m ^Zp{t) 



The right side of this equation is completely independent of the solution of the stationary equation 
we pick. Thus every solution must be a translation of to. i-. 

This concludes our treatment of the Heaviside case. 



5. Convergence for /5 > ^J2\ogK 

In this section, we shall demonstrate convergence for initial data Gq{x) — exp(— e~'^^) with /5 > 
a/2 log K. Our main tool will be Lemma 2.5 along with the knowledge of convergence in the cases 
(3 < ^2\ogK and /3 = oo. 

Let us write Cf for the solution of the recursion relation with initial data Gq{x) — exp(— e"'^^). 
One can check that m^{t) = (G'f)-i(i) is well defined. 
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Lemma 5.1. For /3 > p', G^{x + m^{t)) > Gf{x + m^' {t)) for x > and Cf (x + m^{t)) < 
Cf (x + m^'(t)) for a; < 0. 

Proof: Let us fix to and set Gf (x) = Cf (:c + m'3(to)) and Cf (x) = Cf (x + m^'(to)). Now 
G'J(x) > &^{x) if 

^ /3'm'^'(to)-/3m/'(to) 
X > ^ ^ — - 

and G^{x) < G^\x) if x < ^''»^'faKf"^^(^o) _ Lemma 2.5, there is a Xt so that (x) > (x) 

for X > and G^{x) < Gf{x) for x < Xf Since ^^^(O) = | = G'fJ(O), we see that xt^ — 0. Since 
to was arbitrary, this proves the lemma. □ 

Corollary 5.2. For any e > and large enough t, G^{x + ■m^{t)) < w{x) + e for x > and 
Gf{x + m^{t)) > w{x) — e for X < 0, where w is the limit of the Heaviside case. 

Proof: This follows by setting /3 = cxo in the previous lemma and then using the uniform 
convergence of the Heaviside case.^ 

From now on, we shall write Wc for a traveling wave with speed c and we shall write w for the 
limit in the Heaviside case. 

Lemma 5.3. For ^ < ^/2hgK, Gf (.x + m^(t)) converges to the traveling wave Wc(0) with 
speed c{/3) and normalized to Wc{j3){0) = |. 

Proof: We know that Gf (x + c{/3)t) converges to a traveling wave Wc(/3) with speed c(/3) 
uniformly. Thus Cf (x + m^{t)) = Wc{i3){x + m^{t) — c{j3)t) + o(l). Setting x = 0, we have 
\ = Wc{0){m^{t) — c{/3)t) + o(l). Passing to the limit we see that the limit limt^.oo(m(t) — c{P)t) 
exists. This imphes that 

lim Gf{x + mP{t)) — Wc(8){x + a), 

t— >-oo 

and a is determined by the condition Wc{j3){a) = |. Now of course Wc(/3){x) := Wc{j3){x + a) is a 
traveling wave as well, 



We are now ready to prove convergence for /3 > -\/2 logK. 

Theorem 5.4. For /3 > V21og K, G^{x + m^{t)) converges uniformly to w{x), where w is the 
limit of the Heaviside case. 



Proof: Lemma 5.1 implies that for j3i < < \/2\ogK we have u'c(/3i)(x) < Wc{p^){x) for x > 
and Wc{ii^){x) > Wcds^^^x) for x < 0. Also ^i;c(/3i)(^) — w{x) for x > and Wc(^.)(x) > w{x) for 
X < 0. This implies that we have a pointwise limit 

w{x) — lim Wc{x). 
Since Wc are increasing functions, w is increasing as well. Now we have for small 5 > 



(x - ^2 logic - 5) = y piy)w^/2W+si^ + y)^dy. 



'"'v'Slog K+S 
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Since < < 1, we see that 



^721^+5(3; - ^2 log K) + 0{6) = j p{y)w^/2^+5{x + y)^dy. 

Taking the hmit 5 — )■ 0, we see that 

w{x - \/2\ogK) = j p{y)w{x + y)^dy. 

Since up to translation, this equation has a unique increasing solution and w{0) = |, we see 
that w = w. Now Lemma 5.1 implies that we have the pointwise estimate w^ 2 log K +si^) ~ ^ ^ 
Gf{x + m^{t)) < 'w{x) + e for x > 0, /3 > a/2 log 5 > Q and large enough t. Then taking <5 — !■ 
we see that pointwise G^{x + m^{t)) — )■ w{x) for 2; > 0. We see this in a similar manner for x < 0. 

We can actually extend this pointwise estimate to a uniform one. Let us consider the sequence 
fn{x) = w ^ 2iogA' +j-(^)- know that pointwise /„(x) increases to w{x) for x > and decreases 
to it for X < 0. Also we know that /n(x) — )■ 1 as x — )■ 00 and /n(x) — )• as x — )■ —00. Let e > 0. 
We can take so large that 

< 1 - w{x) < 1 - fn{x) < 1 - /i(x) < e 

for X > A and 

< w{x) < fn{x) < /i(x) < e 

for X < —N. Now < /^(x) < 1 and < /„(x) < 1 for all x. This implies that the sequence (/„) 
is uniformly bounded and equicontinuous. Thus we can pick a subsequence {fnk)k that converges 
uniformly to w on compact sets. We thus have for large enough t and k 

w^/2T^+^{x) - e < Gf (x + m^(t)) < w{x) + e. 

uniformly on [0, 00). Taking A; — > 00 we see that Gt{x + m{t)) — > w{x) uniformly on [0, 00). The 
proof for (—00,0] is similar. ^ 

Knowing convergence to the traveling wave allows us to extract the leading order contribution 
to m'^(t). 

Lemma 5.5. If Gt{x + f{t)) — )■ w{x) uniformly for some f{t), where to is a traveling wave with 
speed c > ■\/2 log K, then /(t + s) — /(t) — )■ cs for any fixed s and — c. 

Proof: First of all we note that a simple indeuction implies that for any two initial data (mea- 
surable and between and 1) Go and Go, sup^ \Gt{x) — Gt(x)| < sup^ |Go(x) — Go(x)|. Also 
Gt{x) = w{x — ct) satisfies f ll.ip . Combining these two remarks, we see that for any fixed s, 
sup^ \ Gt+s{x + fit) + cs) — w{x)\ < sup^ \Gt{x + fit)) — w(x)| — )■ 0. On the other hand, also 
Gt+s(x + f{t + s)) — w{x) uniformly for each s, so comparing the two sequences of functions one 
can argue that f{t + s) — f{t) — )■ cs for each s. This in turn implies that — t- c.^ 

6. The discrete time Brownian bridge and lower order terms for m^(t) 

Up to now, we have showed that Gt{x + m^{t)) converges to a traveling wave uniformly. We 
have also showed that for /3 < ■\/21og~K", m^{t) — c{j3)t — G for some constant G. We know that 
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for /? > log K, the leading order term in m^{t) is \/2\og Kt. So our next goal is to find the 
lower order terms of m^{t) for j3 > \/2\ogK. 

As mentioned in the introduction, Bramson has done this in the continuum time case. In 
discrete time, the problem can be solved with very similar arguments. We shall not repeat all 
of his arguments, but merely formulate the problem in discrete time in a similar manner as the 
continuum problem, prove an estimate that is very important in many other estimates and then 
give a brief sketch of the argument. 

Bramson's main tool in analyzing the KPP-equation is the Feynman-Kac formula. He uses this 

to represent the solution of the equation in terms of an expectation with respect to the Brownian 
bridge. To derive the Feynman-Kac formula in discrete time, let us write Uf = 1 — Gf and iterate 
the recursion relation for U. One obtains 



where ks{y) — logX^jL^^ GsiyY- We can of course interpret Xf — x + Yll=i Vs a random walk 
(actually discrete time Brownian motion). Let us write P"^ for the law of this random walk starting 
at X and for the expectation with respect to it. So we see that 

Ut{x) = (Uo{Xt)e^''=^''*-^^^^^^ . 

We can then split the expectation so that we consider random walks from x to y and average over 
the end point y so we have 

Ut{x)= I ^^e-^t/o(y)E^'^(e^-i'='-(^''))dy, (6.1) 

where y is a random walk from a; to y in t steps with normalized Gaussian increments. is 
the expectation with respect to this random walk and we shall also write P^'^ for the law of it. 

One can check that the density of the joint distribution of (Ifcu ^fej (where < /cj < /cj+i < t) 
is 



where /cq = 0, kg+i = t, yo = x and ys+i = y. Another way to view this process is that it is a 
continuum time Brownian bridge evaluated at integer times. Because of this interpretation, we 
shall write P^'^ for the law of the continuum time Brownian bridge as well and indicate whether 
we are interested in the continuum or discrete time in some other way. 

In analyzing the recursion relation through the Feynman-Kac formula, the main task is to 
determine which paths contribute significantly to the expectation. In the continuum time case, 
Bramson does this by calculating probabilities for the Brownian bridge to hit different sets. This 
involves quite a few technical details and a fair amount of work. In the end, almost all of these 
estimates rely on the fact that one can calculate the following probabihty exactly 



P,^'^(y(s) > for s e [0,t]) = 1 - et. 
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(6.3) 



The proof of this can be found for example in Bramson's work. As mentioned in the introduction, 
the analysis in the discrete time case is formally identical to the continuum time case and we shall 
not be going over the technical details here. What we will do is to demonstrate that one can use 
Fl'^{Y{s) > for s G {0, as in the continuum case. To do this, we shall need the following 

result that Bramson proves. 

Lemma 6.1. Let li, I2 : [0,t] — )■ [—00, 00] be upper semi-continuous functions satisfying li{s) < 
his) for all se[0,t]. Then 



P"'^(y(g) > kjs) for s e [0,t]) 
P"'^(y(s) >/i(s) forse [0,t]) 



is increasing in x and y. 

We note that this result also contains the discrete time case since if / is defined on {0, t}, it can 
be extended to an upper semi-continuous function on [0,t] by setting l{s) = —00 for non-integer 
values of s. 

Using this lemma, we can prove the required estimate in the discrete time case. The approach 
to the proof was suggested by Greg Lawler. For brevity, let us write P^'^(i?o) for P^'^(y(s) > 
for s G {0, in this proof. 

Lemma 6.2. There exists a constant C > so that for x,y > 0, 



P-^(Bo) < C^'+^f + y\ (6.4) 



There is also a constant C" > so that ii x,y > and xy < t, then 



Proof: The lower bound is just an elementary estimate related to the corresponding continuum 
quantity (16.31) : we have 



P^'^(5o) > P"'^(y(s) > for s G [0, t]) = 1 - e 



For the upper bound, we split the random walk into three parts of length |: 
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2 fCXJ fOO 




^0 



X 



2vr(t-2L|j) 



24IJ 



-.e ^L$J r.'^(BQ)dzidz2 
L3J 



1 



2 f'^ 1 (z^ -x)'' 

<V27rte^/ / e'^^F^fUBo) — 

/o Jo L.J ^^(,_2LiJ) 



X 



2-LlJ 



2 



< Ce^^P^ > for s G <^ 0, 



X{s) > for <^ 0, 



By the Gambler's ruin estimate (see tl9j), there is a constant (7 > so that for a: < 

P^(X(s) > for s G {0, < (7^^. 

This gives the desired result for x,y < \/t (the exponential term is bounded for such x and y). 
For x,y > \/t, the upper bound is greater than one so the bound holds in this case as well. Let 
us now consider the < \/i and y > 



By Lemma 6.1, 



Pt'^(F(s) > for s G [0,t]) 



rt'^Y{s) > for s G {0, 
is increasing in x and y. Thus for x < y/t and y > y/i, using (16. 3p we see that 



^xv^r.^ Ft'y(Y(s) > for s G [0, t])) vt, „ ^ 



< 



Pi'^(y(s) > for s G [0,t])) 
l-e-^^(l + x)(l + v^) 



23: 

1 _ e 



C- 



One can then check that for x < \/i and y > \/t, 



2xy 

1 _ ~ 



So we find that 



Pt '^(y(s) > for s G {0, ...,t})<C 



(l + x)(l + i/) 



for all x,y > 0. 
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We shall now briefly go over the final arguments in proving form of the lower order contributions 
to m^{t). We shall completely gloss over the technical details. The case of /3 = \/2 log K and 
(5 > y/2 \ogK need to be treated separately. We shall first consider /3 > ■\/2 log K. 

~ 2^/2 log K 



Lemma 6.3. For /3 > ^/2\ogK, there is a constant C so that m'^{t) > y/2 log Kt——y= — ^ logt+ 



C. 

Proof: The first thing to note is that one can show that there is a constant C so that m'^(t) > 
m^{t) + C", where is the centering term in the f3 = oo case. Thus we only consider the P = oc 
case. We shall also write m{t) = m^{t) for this lemma. The bulk of the technical work in this 
lemma consists of showing that for any fixed yo, x > m{t), y > yo and large enough r (which is 
considered fixed with respect to t), there is a constant C (depending on r) so that 

E^'y (^e^-i^'-(^»)) > C'K*P°'°(F(s) > for s G {r, - r}). 

To do this, one has to work a fair amount to identify the paths with significant weight in the 
Feynman-Kac formula and also show that the measure of this set of paths can be compared with 
X*P°'°(y(s) > for s G {r, ...,t- r}). 

Using (16. 4p . one can then argue that for some constant C (depending on r), 

P°'°(y(s) > for s G {r, ...,t-r})>-. 

i 

These estimates then imply that for any yo, x > m{t), Ut = 1 — Gt, with Gt given by the 
recursion relation with Heaviside initial data, one has 



poo i^^Lll ^ 

Utix) > / Uo{y)^=^K'^dy. 



On the other hand, if one sets x = ^2 log Kt — ■^-^==\ogt + Zi, where l^il < C2\/i for some 
C2 > 0, a quick calculation shows that for any fixed yo < 

Jyo V2nt 

Let us now assume that for any fixed Zi, one can find a t so that if we choose x as above, 
X > m{t). We note that this is equivalent to saying that for any constant C, one can find a t so 
that m{t) < V2 log Kt - ^^^^^^^ logt + C. So for such a t it follows that Ut{x) > de-^'^^^^'^^K 
But with a suitable choice of zi, this will imply that Ut{x) > 1, which is impossible. So we 
conclude that for some constant C, 



m(t) > ^/2log Kt - hgt + C 

2v2 log K 

for all t. As we noted at the beginning, this implies that the same bound holds for m^{t). ^ 
Lemma 6.4. For (3 > a/2 logK, there is a constant C so that ■m'^{t) < -J2 log 2^'i\o^k 1°S^+ 

C. 
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Proof: We consider again U = 1 — G and define Uq{x) = 1 for x < and Uq{x) = Uo{x) for 
X > 0. We note that Uq{x) > Uo{x) for all x and Uq{x) > Uq (x) where Uq is the initial data 
in the Heaviside case, i.e., when f3 = oo. li is then given by the recursion relation with initial 
data Uq, one can check that Uf{x) > Ut{x) and Uf{x) > U^{x) for all t and x. Moreover, Uf is 
strictly decreasing so m*{t) — {U^)~^{\) is well defined. Since is decreasing, m*{t) > m^{t) so 
we only need to show the result for m*{t). Let us write m{t) for the centering term in the j3 = oo 
case and kg for the term in the exponential of the Feynman-Kac formula in the /3 = oo case. 

The majority of the technical work (which we shall skip) for this lemma goes into proving that 

E^'^ (^e^l=ikt-s{Y{s)^ < C"K*Pf""''"(^''^(y(s) > for s e {r, ...,t- r}). 



where x > m(t) + 1, y is arbitrary, C > depending on r, z = x — \/21og Kt + log t + C 
for a suitable C and r is taken large enough (though fixed with respect to t). 

Taking this result as a given, U^{x) > Uf^{x) for all x and t then implies that k^{x) < kt{x) for 
all X and t. Thus we have the same upper bound for the expectation with k* instead of k. Plugging 
this into the Feynman-Kac representation for the recursion relation of and performing some 
estimation, one finds that there is a constant C so that for x > m(t) + 1, 

U;{x) < CK^ / [/o(?/)^=e-^^Pf'^(F(s) > for s e {r, - r})dy. 
Jo \/27rt 

Writing x — zi + \/2 log Kt and using Lemma 6.2 one finds 

lJl{x) < Ce-^^^'' / e^^^^C/o(?/)^=e^^^Pf^(r(s) > for s e {r, t - r})dy 

Jo V27rt 

POO 

<Ct-^l + z)e-^^^'' / e^^^^il + y)Uo{y)dy. 

Jo 

Since /5 > log K, the integral converges and we only care about the part depending on t, 
zi and z. Noting that z — zi — ^^/f^j^logt + C, we see that we can move the f"^ into the 
exponential and get 

U;{x) < C{l + z)e-^^^'. 

Using similar arguments to those in sections 4 and 5, one can show that Ut{x + m*{t)) converges 
uniformly to u = 1 — w. We note that since m*{t) + D > m{t) + 1 for some constant D, we 
can set x = m*{t) + D in this estimate. Since U^{m*{t) + D) converges to u{D) > 0, the 
sequence is bounded from below by some positive number. The inequality above then implies 
that z — m*{t) + D — \/2 log Kt + 2^2\o%k ^'-'S ^ + must be bounded from above (otherwise the 
right side of the inequality could get arbitrarily close to zero) . So we conclude that there is some 
constant C so that 

m\t) < m%t) < y/2log Kt - \ogt + C.^ 



Combining these results shows that for j3 > y/2 logK, 
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m^(t) - a/2 log Kt^ lost 



is bounded. We shall now proceed to the /3 = \/2 log K case. 

Lemma 6.5. Let m : Z+ ^ M satisfy ^ ^2 \ogK as t ^ oo, m{t) > y/2\og Kt - and 
m{t) — m{s) > ^J2 \ogK{t — s) — {t — sf — Ci for sq < s < i and t chosen large enough. Moreover, 
assume that 

ii:*E^(C/o(X(i)); > m{t - s) for s e {0, t - so}) = Dt{z)v{z) 
for some functions v and D^, where z = x — m{t) and Di satisfies 

liminf liminf > 

and So is such that for a fixed m(s) is finite for s e {so,...,t}. Then for /3 = -\/2 logii', 
m(t) — mP{i) is bounded for t > Sq. 

Proof: Consider Mt = max{|m(s) — m^{s)\ : sq < s < t}. The first thing to do in this lemma, 
that would require some work and whose proof we shall skip, is to show that under our assumptions 
concerning m{t), 

E^([/o(X(t)); X{s) > m^jt - s) for s e {0, t - sq}) 
E^([/o(X(i)); > m{t - s) for s e {0, t - Sq}) 

can be estimated from below in terms of a ratio of integrals concerning probabilities of the form 
P"'''(F(t) > for t e {0, s}). Then using Lemma 6.2, one can show that for c > and some 



W{Uo{X{t))- X{s) >mf'{t- s) for s e {0, t - sq}) ^ ^ y{t) 



> C- 



-cMt 



E^{Uo{X{t));X{s) > m{t- s) for s e {0,...,t- So}) " 1 + y{t) + Mfl + z' + Mt 



for some C > if a; > m^{t) + A'^c, where y{t) is some sequence satisfying y{t) — )■ oo as t ^ oo 
and z' = X — m^{t) — Nc- We note that that since Uo{x) > for all x, then always 

K^'E''(Uo(X(t));X(s) > m(t - s) for s e {0, t - Sq}) > 

which implies that v{z) > for large enough z. Indeed taking z large enough and then fixing it, 
our assumptions imply that for large enough t 

K^E%Uo{X{t)y,X{s) >m{t-s) for s e {0, t - Sq}) >C{z), 
where C{z) > 0. Combining these estimates, we see that for large enough (fixed) 2; > and t > 



K'E^Uo{X{t));X{s) > m^{t - s) for s E {0, t - So}) > C^(^) Yip^gy^T^T^ 



e 



Let us now show that m{t) —m^{t) is bounded. We begin by assuming that for some large fixed 
value of t, m^{t) < m{t) and set z > Nc + 1. Then z' > 1 and the right side of the inequahty 
above is at least 
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for a large enough Mt (if no such t can be found, either is bounded or we can find one in the 
case we consider soon). So we conclude that under these assumptions 

K'E''{Uo{X{t));X{s) > m^{t-s) for s e {0, ...,t - Sq}) > C{z)e-'"'^\ 

The second thing that would require some work (which we shall not do) is to show that for 
h < -\/2 log K and large enough (fixed) t (for which m^{t) < m{t) is satisfied), one has 



K^E''{Uo{X{t));X{s) > m^{t-s) for s e {0, t - so}) < C^e 



-bz' 



for X > xi = m^{t) + iV^ for large enough N^. We recall that x — m{t) + z — 171^(1) + Nc + z', so 
z' — z + m{t) — m^lt) — Nc and 

Combining our estimates, we see that 

Now z was large but fixed so if we set c = | and b = 1, then actually for such t, m{t) — m^{t) < 
\Mt + C5. One can than work a bit to show that if m^{t) > m{t), one can perform the same 
arguments with the roles of m(^) and m^{t) switched. 

Combining the estimates for m{t) — m^{t) and m^{t) — m{t), we see that|m(t) — m^(t)| < \Mt + 
Cq. As Mf is an increasing function of t, we see that for all s £ {sq, t}, \m{s)—m^{s)\ < ^M^+Cq 
which implies that Mt < 2Cq. Thus Mt is bounded. ^ 

So all we have to do is to show that m{t) = ^2 log Kt— 2^2iog k ^^S^ + '^t satisfies the conditions 
of the previous lemma with some bounded at- 

Lemma 6.6. For m{t) = ^/2\ogKt - ^^^^^^^ logi + b{t), with 

b{t) = ^^^i-^ log ye^^yUo{y)e-idy^ , 

m^{f) — m{t) is bounded for j3 = \/2 log 

Proof: We shall assume that when talking about b and hence m, we are considering t > 1 so 
that everything is well defined. We begin by noting that b{t) = -^^^= + 0{l). So indeed if wc can 

show that m{t) — m^{t) is bounded we will have showed that m^{t) has the lower order behavior we 
claimed. The form of b{t) also implies that — )■ v^2logX as t — )■ 00 and m{t) > ^2 log Kt — 
for large enough t. Also one can check that m{t) — m{s) > ^/2]ogK{t — s) — {t — sY — Ci for 
s e {so, ■■■,t} and some constants Ci and sq. So we only have to check the main condition of 
Lemma 6.5. 

In this lemma, the majority of work (which we skip) goes into showing that for x > m{t) and 
z — X — m{t) , one has 
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i^*E^([/o(X(t)); X{s) > m(t - s) for s G {0, t - 1}) 

= Ai(z)ir* / f/o(i/)-==e-^P^'^(F(s)>Oforse{0,...,t})d?/, 

jMt V ^TTI 

where Dl{z) — 1 if we first take t — > oo then z — )■ oo. Also — ?■ oo although so slowly that 
it can be replaced by — oo without changing the asymptotic behavior of D^. Due to the simple 
asymptotic form of our initial data and m{t), one can show that the upper limit can be changed 
to t2+'^ for any 5 > without changing the asymptotic behavior of Dj. 

Now for any fixed z > 1 and y E [Mt, t^^^], we can write 

Ft^'{Y{s) > for . e {0, = Fl{z,y)^, 

where by (16 ■4p . {F^)t is bounded from above and below by some positive constants. Also for the 
z, y and t we are interested in, we can write 



V2 log X(z+m(t)) 



where again {F^)t is bounded from above and from below by some positive constants. Plugging 
these into our formula for the expectation, we find 

K*E^(f/o(X(t)); X{s) > m{t - s) for s G {0, t - 1}) 

= Di(z)t-tze-v^^^^e-v'^^("W-v'^^*) [ Uo{y)ye^/^^ye-^dy, 



where Df satisfies liminf2;_!.oo hminff^^oo > 0. Using the definition of m{t), a little more 

estimation then shows that actually 

iC*E^(t/o(X(t)); X{s) > m{t - s) for s G {0, t - 1}) 

poo 2 

= A(2)ri2e-^^^^e-v'^^("W-^^^*) / Uoiy)ye^^^ye-"^dy 
= Dt(z)ze-^^^% 







where liminf^^oo hminfi^oo Dt{z) > and we can use the previous lemma to establish the desired 
result. 

7. Concluding remarks 

We have now concluded our main goal of showing that Gt{x + m^{t)) converges to a traveling 
wave and we have seen the freezing transition at /3 = a/2 log K. Moreover, we have seen that 
in the high temperature case, m^{t) is asymptotically linear in t but at the critical point and at 
lower temperatures there are logarithmic corrections whose coefficients depend on whether or not 
we are precisely at the critical point. 
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In this section, we shall discuss a bit further the relationship of these results to some of the 
problems mentioned in the introduction. First of all, let us see exactly the relationship between 
this problem and that of extreme value statistics (or the maximum of the branching random walk). 

Consider the recursion relation with initial data Gq and let {Xk{t))k denote the branching 
random walk. Consider now the functions 

if* 



\k=l 



We can think of the branching random walk {Xk[t + 1))^ starting at x as consisting of a jump 
from X to X + y and then K branching random walks {Xl{t)) starting at x + y which are all 
mutually independent, identically distributed and depend on y only through the starting point. 
Independence and splitting the expectation into an expectation over y and an expectation over 
the rest of the branching random walks gives 



if' 



Ft+ii^) = I n U.{Xl{t)) dy 

-.Ft{x + y)^dy. 



/27r 

Since Fq = Gq, we see that Ft = Gt for all t. Consider now the /3 = oo case. Then 



if* 



l[Go{Xk{t)) = l{Xk{t) > for all k} 

k=l 

= 1 |mmXfc(t) > o| . 

Since we are dealing with symmetric random variables, we see that this implies that 

Gt{x) = P° (max X kit) < . 

So this is the explicit relationship between our problem in the (3 = oo case and the problem 
of extreme value statistics for certain logarithmically correlated random variables. As discussed 
in [14J, in case the X^ were independent, the limiting distribution is well known and it is of 
the Gumbel form. In fact, the Gumbel distribution is the limit for a large class of dependent 
random variables as well. One specific feature of the Gumbel distribution is that at oo it behaves 
like 1 — Ge~^^. We have seen that for (3 > a/2 log K we have logarithmically correlated random 
variables for which the limiting distribution of the maximum has a tail of the form 1 — Cxe~^^ 
which means they are not in the Gumbel universality class. From the point of view of extreme 
value statistics, the interesting and probably difficult problem is to describe the universality class 
to which our random variables belong. Generalizing our recursion relation might be one way to 
get a hold of some other members in this class. 

The other problem our recursion relation is related to was the problem of the multiplicative 
cascade measures on hypercubes. What we have now showed is that 
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^,(e-^^/3) = G'f(x + m^(t)) = E(exp(-e-^"e-'''"'(*)Z(t))) 

converges as t — )■ oo. Here gt{u, (3) is the Laplace transform of the random variable Mf = 
Z{t). From our analysis, we see that Mf is of the form Att'^e~^'^^^''^Z{t), where a G {0, |, |} 
and At is bounded. Since the Laplace transforms converge to some function which is continuous 
at {gt{u^ (3) — ?• — logw)), we see that Aff converge in distribution to some random variable, 

say . So we have showed that after a deterministic normalization of the multiplicativa cascade 
measures, their total masses converge. 

From Lemma 5.5, m^{t + l)—ml^{t) — )■ c(/3) so this convergence in distribution and decomposing 
Z{t + 1) = e~^^ Tl,k=i Z^if) imply that satisfies the following equation in distribution 

K 

^ g-/9(V+c(/3)) ^ ^/3,fc^ 
k=l 

where 1/ is a standard Gaussian and M^''^ are independent copies of M^^ and also independent of 
V . This equation is of the form that is studied in fT?|. While their main results are about existence 
of solutions to such equations, which is not of direct interest to us since we have a solution, they 
do use a construction that is of interest in our case. Consider the /3 > ^/2\ogK case and let 
g{u, (3) = limj gt{u, (3). Since for (3 > ^/2\ogK, wp = w^ 2 logX ; have for /3 > ^/2\ogK 

log K I \ 

So we see that the Laplace transforms of the limits have a rather simple /3 dependence. A 
similar phenomenon occurs in [12]. Applying their argument in our case implies that is given 
by composing a certain stable process with the random variable M^^^°§^. 

The point of these remarks is to possibly give some tools or approaches to questions such as 
almost sure convergence of Aff or even the existence and structure of the limiting measure. 

Acknowledgments: I wish to thank Antti Kupiainen for introducing me to this subject and 
constantly finding time for useful discussions and advice. I also wish to thank Greg Lawler for 
suggesting the approach to Lemma 6.2 and the referee for helpful comments and suggestions for 
references. Finally I wish to thank the Academy of Finland for financial support. 

References 

[1] L. Addario-Berry and B. Reed. Minima in branching random walks. Ann. Probab., 37(3):1044-1079, 2009. 
[2] I. Benjamini and O. Schramm. KPZ in one dimensional random geometry of multiplicative cascades. Comm. 

Math. Phys., 289(2):653-662, 2009. 
[3] J.D. Biggins. Martingale convergence in the branching random walk. J. Appl. Probability^ 14(l):25-37, 1977. 
[4] E. Bolthausen, J.-D. Deuschel, and O. Zeitouni. Recursions and tightness for the maximum of the discrete, 

two dimensional Gaussian Free Field, http://arxiv.org/pdf/1005.5417v2.pdf, 
[5] M. Bramson. Convergence of solutions of the Kolmogorov equation to travelling waves. Mem. Amer. Math. 

Soc, 44, 1983. 

[6] M. Bramson and O. Zeitouni. Tightness of the recentered maximum of the two-dimensional discrete Gaussian 

Free Field, http://arxiv.org/pdf/1009.3443vl.pdf 
[7] M. Bramson and O. Zeitouni. Tightness for a family of recursion relations. Ann. Probab., 37(2):615-653, 2009. 
[8] D. Garpentier and P. Le Doussal. Glass transition for a particle in a random potntial, front selection in 

nonlinear renormalization group, and entropic phenomena in Liouville and Sinh-Gordon models. Phys. Rev. 

E, 63:026110, 2001. 

24 



[9] A. Champneys et al. Algebra, analysis and probability for a coupled system of reaction-diffusion equations. 
Philos. Trans. Roy. Soc. London Ser. A, 350(1692):69-112, 1995. 
[10] B. Derrida and H. Spohn. Polymers on disordered trees, spin glasses and traveling waves. J. Stat. Phys., 51:817, 
1988. 

[11] B. Duplantier and S. Sheffield. Duality and the Knizhnik-Polyakov-Zamolodchikov relation in Liouville quan- 
tum gravity. Phys. Rev. Lett, 102:150603, 2009. 

[12] R. Durrett and T. Liggett. Fixed points of the smoothing transformation. Z. Wahrsch. Verw. Gebiete, 
64(3):275-301, 1983. 

[13] Y. V. Fyodorov. Multifractality and freezing phenomena in random energy landscapes: an introduction. Physica 
A, 389:4229-4254, 2010. 

[14] Y. V. Fyodorov and J. -P. Bouchaud. Freezing and extreme value statistics in a Random Energy Model with 
logarithmically correlated potential. J. Phys. A, 41(37):372001, 2008. 

[15] Y. V. Fyodorov, P Le Doussal, and A Rosso. Statistical mechanics of logarithmic REM: duality, freezing and 
extreme value statistics of 1/f noises generated by gaussian free fields. J. Stat. Mech., P10005, 2009. 

[16] Y. V. Fyodorov, P Le Doussal, and A Rosso. Freezing transition in decaying burgers turbulence and random 
matrix dualities. Europhysics Letters, 90(60004), 2010. 

[17] Simon C. Harris. Travelling-waves for the FKPP equation via probabilistic arguments. Proc. Roy. Soc. Edin- 
burgh Sect. A, 129(3) :503-517, 1999. 

[18] J. -P. Kahane and J. Peyriere. Sur certaines martingales de Benoit Mandelbrot. Advances in math., 22(2):131- 
145, 1976. 

[19] G. Lawler and V. Limic. Random Walk: A Modern Introduction. Cambdridge University Press, 2010. 

[20] R. Lui. A nonlinear integral operator arising from a model in population genetics, I. monotone initial data. 
SIAM J. Math. Anal, 13(6):913-937, 1982. 

[21] J. Neveu. Multiplicative martingales for spatial branching processes. Seminar on Stochastic processes, Prince- 
ton, NJ:223-242, 1987. 



25 



